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I.  INTRODUCTION 

Research  areas  in  experimental  mechanics  have  been  greatly  expanded  since 
the  introduction  of  the  He-Ne  laser  in  1960.  Initial  engineering  applications 
utilized  holographic  interferometry  to  measure  surface  deformations  of  opaque 
solids.  Original  applications  suggested  that  the  modem  techniques  in  coherent 
optics  would  possibly  change  the  method  of  measurement  of  strains  and  stresses; 
however,  this  has  not  been  the  case,  although  many  engineering  examples  have 
been  successfully  demonstrated.  Some  of  the  authors'  contributions  are  listed 
in  references  [l-7] .  Although  interference  fringes  are  well  defined  in  holog¬ 
raphy,  difficulties  exist  because  of  the  sensitivity  of  the  measurements  which 
require  vibration  isolation  and  the  ability  to  separate  displacement  components 
from  a  single  hologram.  Thus,  holography  has  not  developed  as  a  general  tech¬ 
nique  of  experiemtal  stress  analysis,  and  still  remains  a  technique  of  the 
specialist. 

Objects  illuminated  with  coherent  light  are  observed  to  have  a  granular 
appearance  known  as  the  laser  speckle  effect.  The  details  of  this  structure 
do  not  resemble  the  microscopic  character  of  the  surface,  but  rather  appear  to 
be  random  in  nature.  This  effect  is  described  in  terms  of  the  interference  of 
the  microscopic  nature  of  the  surface  when  illuminated  with  a  laser.  Inter¬ 
ference  of  the  dephased,  but  coherent,  waves  results  in  this  speckle  effect [8]  . 
This  effect  has  provided  for  a  method  of  surface  displacement  measurements 
known  as  laser  speckle  Interferometry,  which  does  not  require  vibration  isola¬ 
tion.  In  addition,  this  technique  is  a  direct  measure  of  in-plane  displacements 
thus,  the  separation  of  components  is  accomplished  easily  [7]  .  All  coherent 
optical  methods  are  a  measure  of  surface  displacements  or  its  derivatives  on 
the  boundary  of  a  solid;  thus,  the  need  still  exists  to  calculate  stresses  and 
strains  at  any  desired  point  in  a  three-dimensional  body  from  the  surface  data. 

The  concept  of  speckle  was  first  applied  to  optical  systems;  however, 
modern  applications  reveal  that  the  speckle  effect  occurs  in  many  areas  of 
physics  and  engineering  [8-10j  as  a  result  of  the  ubiquitous  character  of  the 
random  interference  phenomenon.  Examples  of  analogies  to  the  laser  speckle 
effect  include  radar  astronomy,  synthetic  aperture  radar,  recent  developments 
in  ultrasonic  holography,  and  theory  of  spectral  analysis  of  random  processes. 
Goodman  [8]  illustrates  that  the  basic  concepts  of  speckle  are  much  broader 
than  originally  envisioned  with  the  laser,  and  that  the  results  apply  to  all 
related  phenomena,  so  long  as  the  basic  statistical  assumptions  are  satisfied. 
Because  of  these  parallels,  the  results  of  the  laser  speckle  effects  in  optics 
suggest  the  possibility  of  an  acoustical  speckle  effect  utilizing  ultrasonics. 
The  direct  extension  of  substituting  an  ultrasonic  wave  for  coherent  light  is 
an  obvious  one,  since  it  was  first  performed  in  the  acoustical  holography.  An 
acoustical  speckle  interferometry  study  has  been  demonstrated  to  be  a  direct 
measure  of  object  motion  [l  1,12],  which  closely  parallels  laser  speckle  inter¬ 
ferometry.  The  particular  significance  of  the  ultrasonics  is  the  fact  that 
subsurface  scattering  points  can  be  recorded  and  correlated  to  measure  acous¬ 
tical  signals.  Thus,  acoustical  speckle  interferometry  has  the  potential  to 
measure  displacements  of  interior  points  in  opaque  solids.  The  extension  to 
include  measurements  of  deformable  bodies  will  result  in  a  technique  to  measure 
directly  strains  at  interior  points  in  otherwise  opaque  solids. 
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II .  ELASTODYNAMIC  THEORY 

Acoustical  speckle  interferometry  is  based  locally  on  the  elastodynamic 
response  of  solids  subjected  to  ultrasonic  waves.  The  basic  theory  will  be 
governed  by  the  following  set  of  field  equations  for  a  homogeneous  isotropic 
medium. 

The  stress  tensor  o  (X)  is  related  to  the  gradient  of  displacements  by 
the  following  equation 

"u-‘Vk*U+"(,,i.J+0J,1)  (1> 

where  the  comma  denotes  partial  differentiation  and  X  and  p  are  the  Lame’s 
constants. 

Balance  of  linear  momentum  yields  the  following  equation  of  motion  writ¬ 
ten  in  terms  of  displacement  gradients. 

P  u-  +  (X  +  p)  U  +  pD  =  f  (2) 

i.jj  J.ii  i 

where  p  is  the  mass  density  and  f  are  the  components  of  the  body  force  per 
unit  volume. 

Solution  of  the  equations  of  motion  (2)  must  satisfy  appropriate  boundary 
conditions.  The  displacement  boundary  value  problem  assumes  a  knowledge  of 
displacements  on  the  entire  surface  S. 

U1^)=U1(X)  (3) 

when  X  denotes  the  orthogonal  coordinates  Xj ,  X2 ,  X3  and  Ui  are  known  values 
of  the  surface  displacements.  The  traction  boundary  value  problem  is 

400  -  a  n  =  t±(X)  (4) 

when  the  vector  function  t.(X)  is  prescribed  for  XeS.  Unit  vector  n  is  the 

1  ~  '  j 

outward  normal  vector  for  the  body  P. 

As  an  initial  approximation  to  the  local  response  of  an  elastic  body  to 
an  ultrasonic  wave,  consider  a  plane  displacement  wave  propagating  with  phase 
velocity  C  in  a  direction  described  by  a  propagation  vector  p  as  [13]. 

tj  -  A^ein  (5) 


where 


n  *  K(3>p  -  ct) 


Equation  (5)  describes  a  Diane  harmonic^wave  with  3-?  *  constant  describing  s 
plane  normal  to  the  propagation  vector  p. 

Two  types  of  plane  harmonic  waves  described  by  Equation  (5)  exist  which 
are  solutions  to  Equation  (2) .  Waves  of  the  type  where  3  *  4^  are  called 
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longitudinal  waves  and  the  phase  velocity  is  denoted  as  C  .  Waves  of  this 

-t-  -* 

type,  where  d  •  p  =  0,  are  called  transverse  waves,  and  the  phase  velocity  for 
this  type  of  wave  is  denoted  as  C^. 

Substitution  of  Equation  (5)  into  Equation  (1)  yields  the  following  form 
for  the  components  of  the  stress  tensor. 

°ij  =  [X6ij(d«P£)  +  y(di  Pj  +  djPi)]  ikAeln 
when  the  summation  convention  is  employed. 

The  presence  of  a  subsurface  layer  or  disconunity  generally  produces  a 
change  in  the  wave  propagating  in  a  medium.  Acoustical  speckle  interferometry 
utilizes  the  reflection,  refraction,  and  mode  conversion  of  incident  elastic 
waves  at  the  interface  of  two  elastic  media  as  a  basic  description  of  the 
transmission  and  reflection  of  elastic  waves.  Consider  initially  two  joined 
elastic  half-spaces,  as  illustrated  in  Figure  1.  An  incident  longitudinal 
wave,  denoted  as  a  p-wave,  is  incident  on  the  boundary  of  the  two  elastic 
half-spaces.  At  the  interface,  boundary  conditions  and  conditions  of  conti¬ 
nuity  will  result  in  both  reflection  and  refraction  of  longitudinal  and  trans¬ 
verse  waves  from  an  incident  p-wave.  The  notation  and  description  of  the 
waves  will  be  the  same  as  described  in  reference  [l3]  and  repeated  here  for 
completeness  of  discussion.  Unit  propagation  vectors  are  illustrated  in 
Figure  1,  with  material  properties  of  the  incident  medium  denoted  as  X,  p,  y. 

B  B  B 

Material  properties  of  the  refracted  medium  are  labeled  X,  y  ,  p  .  Super¬ 
scripts,  pn,  on  the  propagation  vectors  are  used  to  describe  the  reflected 
and  refracted  waves  as  given  in  the  following  equations. 

Incident  Longitudinal  wave  (p-wave),  n  -  0. 

-►o 

The  incident  P-wave  will  be  described  by  the  unit  propagation  vector  p  , 

when 


p  =  cos  0  ii  +  sin  0  12 
o  1  o  4 


(7-a) 


Displacement  vector  is  described  as  3°  =  p°;  therefore 


U° 


A  ("cos  0  ii  +  sin  0  t?  I  eiT1° 
o  L  o 


(7-b) 


and 


n  =  K  [Xj  cos0  +  X?  sin  0  -  C  tl 

n1-  n  4  n  n  J 


(7-c) 


Stress  components  for  the  incident  p-wave  are  calculated  from  Equation  (6) 
and  are 


t2?  =  ik  (X  +  2y  cos  0  )  A  eino 

41  O  0  0 


(8-a) 
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o  in 

T22  =  2ik  u  sin  0  cos  0  A  e  o 
o  o  o  o 


(8— b ) 


Reflected  longitudinal  wave,  n  =  1 . 


The  reflected  p-wave  is  described  by  the  unit  propagation  vector  as 


"►l  •>  -V 

p  =  sin  0j  ij  -  cos  0j  i2 


displacement  vector  is  d1  =  p1 


vhere 


U1  =  A}  sin  Ojii-  cos  0ii2 


ni  =  kj  [Xl  sin  ©i  -  X2  cos  0j  -  CjtJ 


Stresses  for  this  wave  are 


121  =  U  sin  20j  Ai  e  i 

t22  =  i  ki  [x  +  2p  cos20}^j  Ai  eiT1l 


(9-a) 


(9-b) 


(9-c) 


(10-a) 


(10-b) 


Reflected  transverse  wave  ( sv-wave) ,  n  =  2. 

Propagation  vector,  displacement  vector  and  stress  components  are 
given  by  the  following  equations. 


p2  =  sin  02ii  -  cos  02i2 


(11-a) 


The  description  of  a  transverse  vertical  wave  is  given  by 
d2  =  i3  +  p2 ,  which  yields  the  following  form  for  the  displacement  vector. 


U2  =  A2  [cos  02"ii~  sin  02"i2Jeiri2 


n2  =  ^2  [xjsin  ©i  -  X2  cos  02-  C2t] 


2  ik  R  y  . 

t21  - - ^ -  cos  2©2A2e  ^2 


t22  =  -ik2ysin  202A2e  2 


(11-b) 


(11-c) 


(12-a) 


(12-b) 
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Refracted  longitudinal  wave,  n  -  3. 

p3  »  sin  03ii  +  cos  03i2 
U3  -  A3  £sin  03ii  +  cos  ©312}  e* 
n 3  ■  K3  |^Xi  sin  03  +  cos  03  -  C3tJ 


m3 


(13-a) 

(13-b) 

(13-c) 


IK  C 

T21  -  — uB  sin  203A3  e' 


(14-a) 


t22  =  ik3^XB  +  2pB  cos2  03J  A3  ein3 
Refracted  transverse  wave,  n  =  4. 


(14-b) 


The  transverse  vertical  wave  in  the  refracted  medium  is  described 


►  O 

by  d3  =  i3Xp  where 


p4  =  sin  04 iy  +  cos  0412 

( 15-a) 

U4  =  A4  cos  04ij  +  sin  04t2Jein4 

(15-b) 

n4  =  K4  Jxj  sin  04  +  X2  cos  04—  C4^.J 

(15-c) 

t2i  =  iK4PB  cos  204A4eiri4 

(16-a) 

4  , 

x22  -  iK4PB  sin  204A4e1Tl4 

( 1 6— b) 

At  the  interface,  a  condition  of  perfect  contact  will  be  assumed  and  the 
displacement  and  traction  vector  at  the  interface  X2  =  0  must  be  continuous 
for  all  Xj  and  t. 


u°  +  u1  +  u2  =  ft3  +  ft4 


t°  +  t*  +  tz  =  t3  +  t 


(17-a) 

(17-b) 


These  equations  written  in  component  notation  are 


U”  +  U!  +  U‘  =  U3  +  it 


(18-a) 
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112  31* 

T2i  +  T2i  +  T2i  =  T2i  +  T2i 


i  =  1,2 


(18-b) 


Equations  (7-c,  9-c,  11-c,  U-c)  require  that  boundary  conditions  are  satislieu 
for  all  Xj  and  t  yields  the  following. 

Kq  sin  Qq  =  Kj  sin  ©i  =  K2  sin©2  =  K3  sin  ©3 


=  K4  sin  ©4 

K0  CL  =  Kj  CL  =  K2  Ct  =  K3  ClB  =  K4  CtB 


(19-a) 


(19-b) 


Equations  (17-a)  and  (7-b,  9-b,  11-b,  13-b)  yield  the  following  conditions  on 
continuity  of  displacements  at  the  interface. 

A0  sin  0g  +  Aj  sin  ©2  +  A2  cos  ©2 

=  A3  sin  03  -  A4  cos  ©4  (20) 

Ao  cos  ©g  -  Ai  cos  ©i  +  A2  sin  02 

=  Aj  cos  ©3  +  A4  sin  ©4  (21) 

Eauation  (18-b),  coupled  with  Equations  (8,  10,  12,  14X  yield  the  following 
conditions  for  conditions  on  the  components  of  the  stress  tensor  at  the 
interface. 


pA0  sin  2©g  -  pAj  sin  2©i  -  cos  2©2 


C  C 

— pB  A3  sin  2©3  -  A4  ~~  yB  cos  2©4  (22) 

CL  Ct 


Kg  +  2u  cos2  0gj  Ag  +  Kj  +  2p  cos2  ©ij  Aj 
-  K2p  sin  2Q2A2  =  K3[xB  +  2pB  cos203]  A3 

g 

+  K4P  sin  204A4 


Equations  (20-23)  can  be  put  in  a  matrix  form  amenable  to  computation 
for  a  system  of  n  layered  media. 


where 


[A]  [X]  =  [B] 


(24) 


M- 


-sinQj 

-COS©2 

sin©  3 

— COS©4 

COS0J 

-sin©2 

COSQ3 

sin©4 

psin2©! 

CL 

V;—  COs2©2 

S 

i^sin2©3 

CL  B 

- „P  COS204 

CI 

(X+2ucos20j ) 

C, 

7^jsin2©2 

S 

ci 

c 

AB+2p 

B  0 

cos^Gz 

CL  B 

— rp  sm2©4 
CT 

_ 

Ai  /Ag 

[x]  = 

A2/Ag 

A3 /Ag 

Ait/Ag 

sin  ©c 

1 

[B]  = 

cos  0g 

p  sin  2©g 

A  +  2pcosQg 

Equation  (24)  forms  the  solution  for  the  amplitude  ratios  Ai/Ag  for  the  reflec¬ 
tion  and  refraction  of  an  incident  P-wave  of  two  perfectly  joined  half-spaces. 
For  a  system  of  n-layers.  Equations  (19)  and  (24)  can  be  utilized  to  determine 
the  propagation  and  energy  partition  at  the  boundary  of  each  layer. 

Some  special  cases  of  the  interfaces  are  of  primary  importance  for 
problems  in  acoustical  speckle  interferometry.  A  liquid-solid  boundary  usually 
forms  the  initial  interface,  and  a  solid-liquid  boundary  forms  the  last  bound¬ 
ary  interface  for  a  solid  immersed  in  a  liquid.  In  addition,  a  subsurface 
layer  from  a  practical  point  forms  a  boundary  of  an  elastic  body  with  air, 
which  is  the  condition  of  a  subsurface  void.  At  this  interface,  the  refrac¬ 
tion  of  elastic  waves  will  be  neglected;  therefore,  only  the  reflection  of 
elastic  waves  at  a  free  surface  will  be  considered.  Since  these  special 
cases  are  of  particular  importance  in  acoustical  interferometry,  they  will  be 
examined  in  some  detail. 
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A.  Liquid-Solid  Boundary 

Suppose  that  the  medium  comprised  of  the  incident  p-wave  is  liquid 
and  the  refracting  medium  is  an  elastic  solid.  For  this  condition,  the 
reflected  SV-wave  is  zero;  therefore,  the  boundary  conditions  reduce  to  the 
following  form. 


U°  +  U1  =  U3  +  U4 
i  i  i  i 

0,1  3.4 

T22  +  T22  =  t22  +  r22 


3,4 
t21  +  T21 


0 


a  =  1,2) 


(25) 

(26) 


Equations  (25)  and  (26)  are  solved  for  the  amplitude  ratios  Ai/A0 ,  and 
are  displayed  in  matrix  form  for  computation. 


— 

r  1 

-  — 

-sin  0 j  sin  63  -cos  84 

Ai  /Aq 

sin  0o 

cos  0j  cos  @3  sin  ©4 

A3/A0 

= 

cos  Oo 

C  r  n  I  C 

-X  XB  +  2v  cos20j  — 5-  p  sin  204 

A4/A0 

A 

CL  L 

_  — 

L_  _ 

B.  Solid-Liquid  Boundary 


The  boundary  between  the  two  layers  is  comprised  of  an^elastic  solid 
reflecting  medium  and  a  liquid  refracting  medium.  The  p-wave  (p°)  represents 
the  incident  wave  and  the  transmitted  sv-wave  is  zero.  Boundary  conditions 
reduce  to  the  following. 


0  12  4 

Ui  +  Ui  +  Ui  =Ui 

(i  =  1,2) 

(28) 

0  1  2  4 
t22  +  t22  +  t22  =  t22 

(29) 

0  ,  1  ,  2 

T2i  +  T21  +  T2l  =  0 

(30) 

Equations  (28,  29)  are  solved,  and  in  matrix  notations  are  written  in 
terms  of  Ai/AO  as 
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-sin  ©i  -cos  02  sin  ©3 

Ai  /Ao 

sin  ©0 

cos  0i  -sin  02  cos  03 

A2/A0 

COS  0q 

C  C 

jx  +  2ucos2©]]  ~psin202  — A^ 

J  T  C 

L 

A3/A0 

A  +  2ucos2©Q 

(31) 


C.  Free  Boundary 

At  a  free  surface,  the  following  boundary  conditions  are  satisfied 
for  the  stress  components. 


0  1  2 
T2i  +  T21  +  T2i  =  0 


i  =  1,2 


(32) 


Solution  of  Equation  (32) 


in  terms  of  the  amplitude  ratios  are  [l3] 


A)  _  sin  20p  sin  20?  -  K2  cos2 20? 
Aq  sin  20q  sin  202  +  K2  cos2 262 


(33-a) 


A]  _  _ 2K  sin  26g  cos  26? 

Aq  sin  20q  sin  2©2  +  K7  cos2202 


(33-b) 


where 

IC 


III.  NUMERICAL  RESULTS 

The  numerical  solution  for  a  layered  elastic  half  space  immersed  in  a 
liquid  medium  involves  the  solution  of  Equations  (24),  (27),  and  (31)  for  the 
amplitude  ratios  Aj /Ao ,  A2/A0,  A3/A0,  and  A*+/A0.  Equation  (27)  is  utilized 
to  determine  the  amplitude  ratios  at  a  liquid-solid  boundary  for  an  incident 
p-wave  in  the  liquid.  At  each  solid  boundary,  Equation  (24)  is  used  to  cal¬ 
culate  the  amplitude  ratios.  At  the  solid-liquid  bound.'  ,  amplitude  ratios 
are  calculated  by  Equation  (31).  Thus,  for  a  multiple-layered  medium,  these 
equations  are  repeatedly  applied  for  the  number  of  interfaces. 

As  an  illustration  of  the  liquid-solid  boundary  equation,  consider  an 
incident  p-wave  in  water  at  a  steel  interface  for  various  angles  of  incidence 
0O .  As  expected,  an  incident  p-wave  gives  rise  to  a  reflected  p-wave  and 
transmitted  p-waves  and  sv-waves  in  the  solid.  Equation  (27)  was  solved 
numerically  for  the  following  data. 


[ 


2(1  -  v) 
1  -  2v 


14 


X  steel 


16.42  x  106  psi 


p  steel 
p  steel 
0  water 
p  water 


0.284  lbm/in3 
11.6  x  106  psi 
0.318  x  106  psi 
0.0361  lbm/in3 


Amplitude  ratios  for  various  values  of  the  incident  angle  ©o  are  given  in 
Figure  2. 

An  incident  p-wave  in  a  solid  at  a  water  interface  gives  rise  to  both 
reflected  p-waves  and  sv- waves  and  a  transmitted  p-wave.  Equation  (31)  is 
used  to  calculate  the  amplitude  ratios  at  this  boundary  and  in  Figure  3  the 
amplitude  ratios  are  shown  versus  the  angle  of  incidence  ©o* 

A  solid  elastic  half  space  for  an  incident  p-wave  in  one  solid  gives  rise 
to  both  reflected  and  refracted  p-waves  and  sv-waves  in  each  solid.  The  amp¬ 
litude  ratios  for  this  boundary  are  calculated  from  Equation  (24).  As  an 
illustration  of  this  solution  procedure  an  incident  p-wave  in  steel  at  a 
steel  aluminum  boundary  is  used  as  data  for  the  calculation  of  the  amplitude 
ratios  for  various  angles  of  incidence  Qq .  Data  for  the  aluminum  are  shown 
in  Figure  4. 


X  Aluminum  =  7.42  x  106  psi 
p  Aluminum  *  4.00  x  106  psi 
p  Aluminum  »  0.100  lbm/in3 

An  example  problem  of  a  two-layered  elastic  half  space  immersed  in  water 
is  chosen  to  illustrate  the  solution  procedure  for  the  layered  medium  problem. 
Geometry  for  this  example  is  shown  in  Figure  5,  and  only  reflected  and  refracted 
p-waves  at  each  interface  are  shown  in  this  illustration;  however,  sv-waves 
exist  at  each  solid  boundary.  The  particular  significance  of  the  p-waves  is 
in  the  application  of  acoustical  speckle  interferometry  to  object  motion 
measurements  in  layered  media.  Reference  [ 12 J  illustrates  that  the  p-waves 
are  amenable  to  object  motion  measurement;  therefore,  this  discussion  will  be 
restricted  to  p-wave  propagation.  Consider  an  incident  p-wave  in  water  at 
interface  1  which  is  a  liquid-solid  boundary.  Solid  1  is  taken  to  be  steel 
with  the  same  properties  for  the  liquid-solid  example.  Superscripts  refer  to 
the  particular  wave  component,  as  discussed  in  the  previous  section,  and  sub¬ 
scripts  refer  to  each  boundary.  A  transmitted  p-wave  at  boupdary3l  will 
become  an  incident  p-wave  at  boundary  2  and  incident  angle  ©2  *  ©i .  For  the 
two- layer  elastic  half  space,  the  incident  waves  at  each  boundary  are  listed 
in  the  following: 


Boundary  1: 

?? 

=  P° 

,  ©1 

-  0° 

Boundary  2: 

-*•0 

P2 

-*■3 

“  Pi 

,  ©$ 

“  9? 
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AMPLITUDE  RATIO 


ANGLE  OF  INCIDENCE  eQ  (DEGREES) 


Figure  2.  Amplitude  ratio  of  a  water  steel  interface 
for  an  incident  p-wave  in  water. 
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AMPLITUDE  RATIO 


R 


Figure  3.  Amplitude  ratio  of  a  steel  water  interface 
for  an  incident  p-wave  in  steel. 
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AMPLITUDE  RATIO 


0  2  4  6  8  10  12  14  16  18  20  22 


ANGLE  OF  INCIDENCE  GQ  (DEGREES) 


Figure  4.  Amplitude  ratio  of  an  aluminum  steel  elastic  half  space 
for  an  incident  p-wave  in  aluminum. 
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I 


Figure  5.  Reflection  and  refraction  of  p-waves  in  a  layered  medium 
for  an  incident  p-wave  in  the  liquid  medium. 
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Boundary  3: 

-+0 

P3 

■+1 
=  P2 

,  ©3 

=  0° 

Boundary  4: 

-vO 

p4 

*  p2 

0 

,  04 

=  02 

Boundary  5: 

-*■5 

p0 

- 

,  e> 

“  04 

Boundary  6: 

-*0 

p6 

-  Ps 

. 

=  ©5 

Amplitude  ratios  at  each  interface  for  an  initial  incidence  angle  0°  =  0° 
are  listed  in  the  following  for  the  two-layer  elastic  half  space.  The  com¬ 
puter  program  for  these  calculations  is  listed  in  Appendix  A. 


Liquid-Solid  1  Boundary  1 
Aj/Aq  =  0.93810 

A3/A0  =  0.061905 

A4/Ao  =  0.00 


Solid  1 -Sol id  2  Boundary  2 


A^Aq  =  -0.45965 

A2/Aq  =  0.00 

A3/Aq  *  1.4597 

A^/Aq  =  0.00 


Solid  1-Liquid  Boundary  3 
A^Aq  *  -0.93809 

A2/Ag  =  0.00 

A3/Aq  -  1.9381 


Solid  2-Liquid  Boundary  4 
Aj/Aj,  =  -0.8411 

A2  /  Aq  —  0 . 00 

A3/Aq  =  1.8411 


Solid  2-Solid  1  Boundary  5 


Aj/Aq  =  0.45965 

A2 / A^  =  0 . 00 

Aj/Aq  =  0.54035 

A^/Aq  -  0.00 


Solid  1-Liquid  Boundary  6 
Aj/Ac  -  -0.93809 

A2 /Aq  =  0 • 00 

Aj /Aq  =  1.9381 
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Amplitude  ratios  at  each  interface  for  an  initial  incidence  angle 
0°  *  2Cr  are  listed  as  follows: 


Liquid-Solid  Boundary  1 


Ai  /Ao 

=  0.92684 

A3/A0 

=  0.17302 

A4/A0 

=  0.65243 

Solid  1-Solid  1  Boundary  2 

Ai  /Aq 

=  0.44433 

A2/A0 

=  0.09316 

A3  /Ao 

=  1.4611 

A4/A0 

=  0.017700 

Solid  1-  Liquid  Boundary  3 

Ai  /Ao 

=  0.81650 

A2/Aq 

=  0.69193 

A3/A0 

=  2.0430 

Solid  2-  Liquid  Boundary  4 

Ai  /Aq 

=  0.72645 

A2/Aq 

=  0.65056 

A3/A0 

=  1.9423 

Solid  2-  Solid  1  Boundary  5 

Ai  /  Ao 

-  0.44651 

A2/Aq 

=  0.08059 

A3 /Aq 

=  0.53946 

A4/A0 

■  0.0077461 

Solid  1-  Liquid  Boundary  6 

Ai/Aq  *  -0.81650 
A2/Ao  =  0.69193 
A3/A0  -  0.02430 

IV.  ACOUSTICAL  SPECKLE  INTERFEROMETRY 

As  a  tool  of  quantitative  nondestructive  testing  of  materials,  the 
spectral  analysis  of  ultrasonic  pulses  in  elastic  solids  has  attracted  wide 
attention  in  recent  years.  Krautkr&ner  [14]  suggested  in  1959  that  measure¬ 
ment  of  the  pulse  shape  could  be  used  to  determine  the  size  and  orientation 
of  flaws  in  materials.  Methods  for  determining  the  spectral  content  of 
broad-band  ultrasonic  pulses  reflected  and  scattered  in  materials  were  subse¬ 
quently  developed  by  Gericke  [l5].  In  his  pioneering  work,  he  showed  that 
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the  spectral  content  of  pulses  propagating  in  steel  specimens  was  greatly 
dependent  on  the  material's  microstructure.  In  later  work  [l6],  he  investi¬ 
gated  the  spectral  content  of  echoes  from  cylindrical  cavities  oriented 
perpendicular  and  parallel  to  the  direction  of  wave  propagation.  From  these 
experiments  he  concluded  that  it  might  be  feasible  to  use  ultrasonic-spec¬ 
troscopy  measurements  in  conjunction  with  normal  pulse-echo  techniques  to 
determine  not  only  the  location  of  a  defect,  but  also  its  shape.  Brown  [17] 
has  reviewed  some  results  obtained  by  Lloyd  [l8]  on  the  spectral  response  of 
various  conical  and  flat-topped  cylinders,  which  could  be  explained  quali¬ 
tatively  in  terms  of  the  cylinder's  dimensions  and  geometry.  Sachse  [l9] 
has  reported  an  experimental  investigation  of  the  spectral  analysis  of  wide¬ 
band  ultrasonic  pulses  as  scattered  by  a  circular,  cylindrical,  fluid-filled 
cavity  in  an  elastic  solid.  It  was  shown  that  either  the  diameter  of  the 
cavity  or  the  wave  speed  of  the  fluid  can  be  determined  from  the  time  record 
of  the  scattered  pulses.  Furthermore,  an  empirical  formula  was  developed  to 
relate  the  diameter  and  the  wave  speed  to  the  intervals  between  the  minima 
or  maxima  of  the  power  spectra.  The  author  concluded  that  the  method  of 
ultrasonic  spectroscopy  which  are  originated  by  Gerick  can  be  used  to  detect 
the  dimensions  and  mechanical  properties  of  flaws  and  impurities  inside 
elastic  solids. 

Theories  of  the  scattering  of  acoustic  waves  have  been  extensively 
devleoped  since  the  fundamental  work  of  Lord  Rayleigh  [20] .  A  review  of  the 
existing  theories  as  well  as  some  numerical  results  of  scattered  waves  in 

elastic  solids  have  been  published  in  the  monograph  by  Pao  and  Mow  [2l].  The 
case  of  a  circular,  cylindrical  cavity  inside  an  elastic  solid  was  investi¬ 
gated  in  detail  by  White  [22],  He  measured  amplitudes  and  arrival  times  of 
plane  harmonic  waves  with  long  duration,  using  standard  ultrasonic  techniques. 
His  measured  angular  distributions  of  the  scattered  wave  amplitudes  for  both 
incident  longitudinal  and  shear  waves  agreed  very  well  with  his  theoretical 
analysis. 

The  interaction  of  sound  waves  and  pulses  with  solid  cylinders  immersed 
in  a  fluid  has  been  of  interest  in  underwater  acoustics  for  some  time.  Most 
extensively  studied  has  been  the  scattering  of  pulses  incident  normal  to  the 
axis  of  the  cylinder.  In  1967,  Goodman,  Bunney,  and  Marshall  [23]  reported 
the  existence  of  a  Rayleigh-type  circumferential  wave  propagating  around  a 
solid  cylinder.  Later,  Neubauer  [24-26]  reported  the  results  of  experiments  in 
which  pulses  were  used  to  study  the  propagation  of  the  circumferential  waves  on 
aluminum  cylinders  immersed  in  water.  In  1969,  Bunney,  Goodman,  and  Marshall 
[27]  also  reported  observation  of  the  Rayleigh-type  wave.  Several  experiments 
of  the  scattering  of  broad-band  ultrasonic  pulses  by  cylindrical  inclusions  in 
elastic  solids  have  been  performed  by  Sachse  and  Chian  [28]  and  Bifulco  and 
Sachse  [29].  The  inclusions  were  fluids  or  solids  whose  wave  speed  was  less 
than  that  of  the  solid  matrix  material.  It  was  shown  that  in  a  non-destructive 
testing  application,  the  size  and  wave  speed  of  the  inclusion  could  be  deter¬ 
mined  from  measurements  of  the  arrival  times  of  various  scattered  pulses.  In 
a  recent  work  [30]  measurements  are  reported  of  the  arrival  times  of  broad-band 
ultrasonic  pulses  scattered  by  a  circular,  cylindrical  solid  inclusion  imbedded 
in  a  matrix  whose  longitudinal  wave  speed  is  lower  than  that  of  the  scatter. 

In  many  instances  in  practice,  it  is  highly  desirable  to  obtain  an 
accurate  and  rapid  measure  of  the  stress  in  a  structural  member.  Ultrasonics 
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offers  a  great  tool  for  the  realization  of  a  system  capable  of  performing  this 
task.  Ultrasonic-pulse  spectroscopy  measurements  has  been  used  to  measure  the 
stress-induced  interference  effects  between  two  shear  waves  propagating  in 
uniaxially  deformed  specimens  of  aluminum  [3l].  This  technique  shows  consid¬ 
erable  promise  as  a  means  for  measuring  and  monitoring  the  applied  stresses 
in  materials. 

Ensminger  [32]  documents  many  modern  applications  dealing  with  topics 
such  as  thickness  measurements,  inspection  of  metals  and  non-metals  and  the 
determination  of  bond  integrity.  Ultrasonic  inspection  is  further  being 
extended  to  tissue  examinations  and  bone  scans  in  the  medical  field  [33]. 

As  an  experimental  nor  -destructive  stress  analysis,  the  polarized  trans¬ 
verse  ultrasonic  waves  have  provided  a  new  method  of  measurements  which  is 
christened  the  acoustoelasticity.  Similar  to  the  photoelastic  effect,  the 
velocity  of  the  transverse  soundwaves  traversing  a  stressed  elastic  solid  is 
not  a  constant  as  it  would  be  in  an  isotropic  solid.  The  velocity  depends  on 
the  direction  of  the  particle  motion  (polarization),  the  direction  of  wave 
propagation,  and  the  state  of  stress.  With  the  recent  advancement  of  ultra¬ 
sonic  techniques,  one  can  generate,  transmit  and  detect  plane  poLarized  ultra¬ 
sounds.  Therefore,  stress-induced  acoustical  birefringence  can  be  measured. 

As  in  photoelasticity,  the  difference  of  propagation  velocities  in  shear  waves 
polarized  in  principal  stress  directions  is  proportional  to  the  principal- 
stress  differences  in  the  plane-stress  state.  Thus,  this  technique  is  a 
nondestructive  test  method  and  determines  stresses  by  ultrasonic  shear  waves 
which  do  not  require  a  transparent  plastic  model  [34-36]. 

Plastic  deformation  in  crystalline  solids  by  means  of  high-f requepcy 
stress  wave  propagations  has  been  studied  by  a  number  of  investigators  [37,38]. 
The  study  of  the  relation  between  the  propagation  of  ultrasonic  waves  and  plas¬ 
tic  deformation  in  single  crystals  has  been  discussed  in  detail.  Specifically, 
changes  in  attenuation  and  velocity  of  ultrasonic  waves  have  been  measured  as 
a  function  of  plastic  strain  during  the  deformation  process.  These  techniques 
have  also  been  applied  to  the  investigation  of  some  aspects  of  strain  hardening 
mechanisms  [39]. 

Knowledge  and  control  of  the  transmission  and  reflection  of  ultrasonic 
waves  in  layered  media  is  an  important  basic  problem  in  underwater  acoustics. 
The  problem  of  transmission  of  sound  at  oblique  incident  angles  through  a 
solid  multilayer  system  of  plane  parallel  plates  was  treated  by  Thompson  [40], 
Brekhovskikh  [41],  and  others.  The  most  detailed  solution  of  the  multilayer 
problem  is  by  Shaw  and  Bugl  [42],  whose  theoretical  approach  was  directed 
toward  the  investigation  of  conditions  under  which  interface  waves  may  be  gen¬ 
erated.  In  a  recent  work  [43],  the  wave  equation  was  solved  to  determine  the 
transmission  and  reflection  coefficients  for  plane  waves  at  oblique  incidence 
on  a  system  of  n  layers  of  plane  parallel  plates. 


Acoustic  holography  is  the  newest  branch  in  acousting  imaging  technology. 
It  is  directly  comparable  with  optical  holography  in  the  phenomenological 
sense,  because  to  make  a  hologram  it  is  necessary  that  the  two  beams  of  radia¬ 
tion  used  have  the  ability  to  interfere  with  one  another.  Acoustic  holography 
underwent  very  rapid  development  during  the  last  few  years,  and  simulated  the 
general  field  of  acoustic  imaging  through  the  introduction  of  novel  techniques 
for  visualizing  acoustic  excitation.  In  contrast  to  optical  holography,  which 


has  sometimes  been  described  as  a  solution  looking  for  a  problem,  acoustical 
holography  was  from  the  outset  developed  for  a  specific  purpose.  This  purpose 
is  that  of  all  acoustical  probing  systems,  namely,  detecting,  locating,  and 
imaging  a  structure  immersed  in  a  medium  opaque  to  electromagnetic  radiation. 
This  description  encompasses  a  variety  of  specific  problems,  such  as  non¬ 
destructive  testing,  medical  imaging,  oil  exploration,  and  underwater  imaging. 

A  recent  development  in  laser  technology,  referred  to  as  speckle  inter¬ 
ferometry,  has  provided  for  a  new  technique  in  experimental  mechanics  to 
measure  surface  displacements  of  deformed  solids.  Basically  this  method 
records  the  laser  intensity  variation  across  the  surface  of  an  object  due  to 
the  random  interference  of  the  microscopic  character  of  the  illuminated  sur¬ 
face.  Figure  6  illustrates  the  basic  method  for  recording  a  double  exposure 
laser  speckle  photograph.  A  high  resolution  photographic  film  records  the 
granular  pattern  on  the  film  plane  of  a  camera.  If  the  object  is  displaced 
from  a  reference  configuration,  the  speckles  are  mapped  to  a  new  location  on 
the  film  plane  corresponding  to  the  object  motion.  Data  analysis  consists 
of  pointwise  filtering  as  illustrated  schematically  in  Figure  7.  This  method 
of  data  analysis  assumes  a  constant  value  of  displacement  within  the  small 
area  of  illumination,  and  the  resulting  interference  pattern  is  obtained  by 
taking  optically  the  Fourier  transform  of  the  amplitude  transmission  function 
of  the  photographic  transparency.  Interference  fringes  are  therefore  a  direct 
measure  of  local  object  motion  within  the  area  of  illuminat ion .  The  discrete 
global  response  of  the  body  is  obtained  by  pointwire  filtering  of  many  small 
areas  of  the  photograph  transparency. 

Because  the  speckle  effect  is  statistical  in  nature  and  the  basic  con¬ 
cepts  occur  in  many  areas  of  physics  and  engineering,  the  previously  devel^^d 
techniques  in  optics  suggest  a  parallel  development  in  ultrasonics.  In  a 
recent  work  [11,12]  a  basic  theory  of  one-dimensional  pulse-echo  and  continu¬ 
ous  wave  acoustical  speckle  interferometry  was  developed.  This  method  utilized 
the  statistical  concepts  of  speckle  and  one-dimensional  correlation  of  refer¬ 
ence  and  displaced  signals  to  determine  object  motion  of  subsurface  scattering 
layers.  This  study  established  the  feasibility  of  the  acoustical  speckle 
concept  for  object  motion  measurements. 

V.  BASIC  THEORY  OF  DISPLACEMENT  MEASUREMENTS  UTILIZING  ULTRASONICS 

This  section  presents  the  geometrical  concepts  of  displacement  functions 
of  a  continuum  and  the  correlation  of  a  reference  and  deformed  ultrasonic 
signal  as  a  measure  of  object  motion.  The  basic  problem  may  be  described  in 
the  following  way:  given  the  measured  ultrasonic  global  response  of  signals 
before  and  after  deformation,  it  is  then  required  to  calculate  the  local 
object  displacement  vector  from  the  global  reference  and  deformed  configurations 
of  a  body.  Thus  the  basic  problem  is  geometrical  in  nature  and  initially  is 
not  concerned  with  the  external  forces  necessary  to  cause  the  geometrical 
deformation. 

Let  the  position  of  the  points  in  a  continuum  in  the  initial  (unstrained 
state)  be  given  with  respect  ot  the  coordinates  Xi ,  X2 ,  X3 ,  as  illustrated  in 
Figure  8.  The  body  is  deformed  with  respect  to  the  unstrained  coordinates 
and  points  are  displaced  to  new  coordinates  denoted  as  Xj ,  X2,  X3  (Figure  8) . 


Figure  6.  Formation  of  a  laser  speckle  photograph. 
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Data  analysis  for  pointwise  filtering  in  in-plane  measurements. 


The  coordinates  of  an  arbitrary  point  in  the  deformed  configuration  are 
expressed  in  terms  of  the  original  coordinates  by  the  following  equation. 


h  =  Xi  +  Ui(Xl •  x2*  x3.  t) 


(34) 


where  are  called  displacement  components  and  are  taken  to  be  continuous 

functions  of  the  original  coordinates  Xj ,  X2 ,  X3 .  Thus,  in  general,  the 
displacement  functions  vary  from  point  to  point  in  a  deformed  body. 

Suppose  a  continuous  mode  or  pulse  echo  ultrasonic  scan  is  made  of  a 
rough  scattering  surface  with  reference  area  M  as  illustrated  in  Figure  9. 

The  field  reflected  from  the  scattering  surface  at  the  plane  of  the  receiving 
transducer  in  a  plane  immediately  adjacent  to  that  surface  is  described  by  a 
complex  function  a(S,n),  which  represents  a  reflected  p-wave  or  SV-wave  from 
an  incident  p-wave  in  an  elastic  solid.  The  complex  field  P(X)  where  X  =  (X},X2) 
denotes  the  ordered  pair  of  numbers  in  the  plane  of  the  receiving  transducer 
which  is  parallel  to  (C,n)  and  represents  the  amplitude  of  the  reference 
signal,  which  is  the  field  of  interest  from  the  acoustical  signal.  Details 
of  both  pulse  echo  and  continuous  mode  scanning  techniques  of  scattering  sur¬ 
faces  are  described  in  reference  [ll,12]. 

Now  suppose  that  the  receiving  transducer  is  returned  to  the  original 
position  and  the  scattering  surface  is  displaced  from  its  reference  configura¬ 
tion.  A  complex  field  P(X')  represents  the  amplitude  of  the  received  signal 
where  the  coordinates  X"  *  Xi',X2  are  parallel  to  the  plane  (£,n).  If  the 
displacements  functions  in  Equation  (34)  are  substituted  into  the  representa¬ 
tion  of  the  displaced  signal,  the  special  form  reduces  to  the  following: 

pqr)  *  p(Xj  +  Uj ,  x2  +  u2)  (35) 


where  the  received  signal  is  assumed  to  be  a  function  of  the  coordinates 
parallel  to  the  plane  of  the  scattering  surface.  The  reference  and  displaced 
signals  are  illustrated  in  Figure  10. 

The  problem  now  becomes  one  of  correlating  the  signals  P(X)  and  P(x') 
to  measure  the  local  displacement  components  Uj  and  U2.  The  basis  for  this 
measurement  is  the  autocorrelation  function  of  the  two  signals  which  is  given 
in  the  following  equation: 

C(U!,U2)  *  J"  P(Xj ,X2)  P(Xi  +  Ui ,X2  +  U2)  dXidX2  (36) 

M 

The  one-dimensional  correlation  of  continuous  mode  scanning  for  a  rigid 
deformation  is  discussed  in  reference  [12].  This  study  demonstrated  that 
Equation  (36)  forms  the  basic  concept  for  object  motion  measurement  of  sub¬ 
surface  acoustic  scattering  layers  which  agrees  with  the  statistical  rela¬ 
tionships  in  laser  speckle  interferometry. 

Although  Equation  (36)  has  been  demonstrated  to  be  a  measure  of  object 
motion,  some  limitations  are  imposed  on  the  applications  because  the 
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Figure  9.  Free  space  geometry  for  acoustical  speckle. 
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Figure  10.  Reference  and  deformed  surface  for  correlation 
of  acoustical  signals. 
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displacement  functions  Uj  and  U2  are,  in  general,  continuous  functions  of  the 
initial  coordinates  (Xi ,X2)  and  are  the  unknown  quantities  to  be  determined. 

If  the  displaced  surface  is  just  a  uniform  translation  in  the  Xj ,X2  directions, 
then  Ui  ,U2  are  constant  and  the  application  of  Equation  (36)  is  straightforward; 
however,  in  the  presence  of  a  deformed  surface,  Uj  and  U2  cannot  be  interpreted 
in  terms  of  a  global  uniform  translation.  The  basic  approach  used  in  computing 
the  correlation  of  P (3C)  and  P(X')  is  to  assume  the  basic  theory  of  data  analysis 
in  pointwise  filtering  used  in  single  beam  laser  speckle  interferometry 
(Figure  11).  Object  motion  displacements  are  obtained  by  passing  a  laser  beam 
through  a  small  area  of  a  double  exposed  transparency  which  contains  both  the 
reference  and  deformed  images.  Interference  fringes  are  observed  as  intensity 
of  the  light  amplitude  in  the  transform  plane.  This  analysis  utilizes  the 
shift  theorem  which  restricts  the  displacement  components  to  be  constant  with¬ 
in  the  area  of  illumination.  Practical  applications  for  a  He-Ne  laser  restricts 
the  area  to  be  a  circular  region  approximately  2  ram  in  diameter. 

This  restriction  in  optics  suggests  the  following  approach  in  acoustical 
speckle  interferometry.  Let  the  reference  scan  be  recorded  over  some  reference 
area  M  in  the  Xj  ,X2  plane,  as  illustrated  in  Figure  11.  Suppose  that  within 
seme  small  neighborhood  of  a  point  eP(X  )  the  displacement  is  assumed  to 

be  uniform.  The  acoustical  signal  corresponds  to  a  small  neighborhood  of  the 
point  Pq  eP(X).  Therefore  the  Equation  (36)  can  now  be  used  in  the  form 

AC(U!,U2)  =  J  AP(X2 ,X2)AP(X!  +  Uj  ,X2  +  U2)  dXjdX2  (37) 

M 

where  AP(X)  and  AP(X") correspond  to  a  small  surface  surrounding  the  points  of 
interest.  The  correct  values  of  Uj  ,U2  corresponding  to  the  object  motion 
displacement  will  result  in  a  maximum  value  of  the  correlation  function 
AC(Uj ,U2).  This  correlation  procedure  corresponds  to  the  data  analysis  at 
pointwise  filtering  in  optics  and  thus  completes  the  analysis  of  the  displace¬ 
ment  of  a  point.  The  procedure  to  determine  the  discrete  displacement  values 
will  follow  the  procedure  in  optics  as  discussed  in  reference  [7]. 

A  discussion  of  the  cross-correlation  peak  should  include  the  effects 
of  deformation  on  the  local  accuracy  of  the  correlation  function.  This  prob¬ 
lem  in  local  deformation  of  an  object  is  similar  to  the  problem  of  image 
correlation  with  geometric  distortion  on  the  accuracy  of  image  restoration 
[44].  Basically,  the  image  restoration  is  investigated  for  distortions 
represented  by  an  affine  transformation  of  image  coordinates,  which  is  an 
offset  plus  scaling  and  rotation  of  each  image  coordinate  axis  [45-47].  The 
image  analysis  assumes  a  fixed  geometric  distortion  and  then  a  proper  choice 
of  integration  area  in  the  cross-correlation  can  be  optimized  for  minimum 
local  error.  Geometric  distortion  for  deformable  bodies  is  generally  unknown; 
therefore,  this  represents  the  inverse  problem  to  image  registration.  A 
discussion  of  this  problem  will  include  the  restrictions  and  assumptions  in 
the  linear  theory  of  elasticity  and  an  example  problem  will  be  chosen  to 
illustrate  the  effects. 

Let  a  point  P  be  displaced  to  a  point  ?' ,  and  point  Q  displaced  to  Q'. 

The  displacement  of  point  Q  can  be  expanded  in  a  linear  Taylor's  series  as 
expressed  in  the  following  form. 
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Reference  and  displaced  subsets  of  the  original 
and  displaced  surfaces. 
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Equation  (38)  can  be  expressed  in  terms  of  the  infinitesimal  strains  and 
rotation  components  which  results  in  the  following  expression  £483. 

Uj  |  Q  =  Uj  |  P  +  AUX  (39) 


U2  [  Q  =  U2  |  P  +  AU 2 


where 


AUj  =  en  AXj  +  (ei2  -  w3) 


AU2  =  (ei2  +  io3)AXl  +  £22^x2 


The  correlation  function  can  now  be. put  in  the  following  form. 


AC(Ui  +  AUi,  U2  +  AUZ)  =  JaP(X!,X2)  AP(Xj+  Uj  +  AUj ,X2 


+  AU2)  dx  dx^  (40) 

Equation  (40)  has  been  demonstrated  to  be  a  measure  of  object  motion, 
taking  into  consideration  the  restrictions  and  assumptions  in  the  linear  theory 
of  elasticity;  however.  Equation  (36)  represents  the  general  measure  of  object 
motion.  In  the  presence  of  a  largely  deformed  surface  U3  and  U2  can  be  inter¬ 
preted  in  terms  of  a  rigid  body  translation  and  rotation  and  normal  and  shear 
strains.  Several  numerical  examples  will  be  discussed  in  order  to  illustrate 
the  application  of  numerical  correlation  in  object  motion  mepsurements. 
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VI.  NUMERICAL  EXAMPLES  AND  CONCLUSIONS 

Several  numerical  examples  will  be  discussed  to  illustrate  the  application 
of  the  numerical  correlation  in  object  motion  measurements.  A  reference  signal 
P(X)  corresponding  to  the  undeforaed  scan  is  illustrated  in  Figure  12.  A  line 
along  the  Xj  axis  is  shown  in  Figure  13.  The  reference  area  M  for  all  exam¬ 
ples  presented  will  comprise  an  array  of  100  x  100  data  points.  The  displaced 
signal  P(X")  will  be  correlated  with  the  reference  signal  for  the  specific  cases 
of  displacement  components  Uj  and  U2 .  A  correlation  computer  program  is  listed 
in  Appendix  B. 

For  a  case  of  uniform  translation  in  the  absence  of  geometrical  distortion, 
some  simplifications  in  the  integration  area  M  allowed  for  a  considerable  reduc¬ 
tion  in  computer  time.  The  procedure  utilized  was  to  perform  the  correlation 
over  two  orthogonal  lines  of  data  points  which  comprise  a  length  where  object 
motion  is  assumed  to  be  uniform.  In  the  examples  discussed,  the  five  data  points 
for  this  interval  were  used.  Figure  14  is  a  graphical  display  of  three  examples 
of  uniform  translation.  Figure  15  shows  the  surface  P(X')  for  a  uniform  trans¬ 
lation  in  the  X* ,X2  directions.  In  order  to  demonstrate  the  generality  of  the 
computer  program,  an  example  of  non-uniform  translation  without  geometric  dis¬ 
tortion  is  discussed.  The  displacement  along  a  line  y  =  sin  x  for  0  £  *  tr 
is  chosen  to  illustrate  non-uniform  translation  and  the  results  are  shown  in 
Figure  16. 

Several  examples  for  the  case  of  uniform  and  non-uniform  translations  in 
the  presence  of  geometric  distortion  (deformation)  were  solved.  The  results 
illustrated  that  for  infinitesimal  deformation  (linear  elasticity),  geometric 
distortion  had  no  effect  on  the  correlation  between  reference  and  deformed 
signal . 

Several  examples  of  large  deformations  on  the  correlation  of  reference 
and  deformed  signals  were  considered  to  determine  the  accuracy  of  Equation  (40) 
for  large  strains.  Figures  17  through  19  are  examples  of  the  deformed  surface 
P(X')  corresponding  to  large  strains  without  rigid  body  motion.  A  comparison 
of  these  surfaces  with  the  reference  surface,  as  illustrated  in  Figure  12, 
shows  the  effect  of  geometric  distortion  relative  to  the  reference  surface. 

In  each  example  of  large  geometric  distortion  the  correlation  accuracy  was  the 
same  as  the  cases  considered  for  uniform  translation  and  infinitesimal  defor¬ 
mations. 
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re  13.  Reference  signal 


Figure  14.  Uniform  translation  of  a  point  P  initially- 
located  at  the  origin. 
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Figure  15.  Acoustical  surface  corresponding  to  a  uniform  translation 
without  geometric  distortion  -5  units  Xj  direction 
-8  units  in  X2  direction. 


Figure  17.  Acoustical  surface  P(X')  corresponding  to  a 

large  strain  ev  =  0.2. 

X1 


Figure  18.  Acoustical  surface  P(X')  corresponding  to 

large  strain  e  =  0.2. 

x2 


Figure  19.  Acoustical  surface  P(X ')  corresponding  to  a 
large  strain  y  =  0.2. 
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APPENDIX  A 

PROGRAM  TO  CALCULATE  THE  AMPLITUDE  RATIO 
FOR  A  TWO-LAYER  ELASTIC  HALF  SPACE  IN  A  LIQUID 
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PROGRAM  TO  CALCULATE  THE  AMPLITUDE  RATIO 
FOR  A  TWO  LAYER  ELASTIC  HALF  SPACE  IN  A 
LIQUID 

EACH  NUMBER  CORRESPONDS  TO  AN  INTERFACE 

PROPERTIES  OF  THE  LIQUID 
XLAL=318000 . 

RHL= . 0361 

PROPERTIES  OF  SOLID  1 
XLA1=16420000. 

XM1=1 1600000 . 

RH1= . 284 

PROPERTIES  OF  SOLID  2 
XLA2=7420000. 

XM2=4000000 . 

RH2=. 100 
CLL=4S63 • 

CL2=20341 . 

CT2=10360 . 

CL1=19348. 

CT1=10469. 

LAYER  1 - LIQUID-SOLID  1  INTERFACE 

THETA0=20, *3. 1416/180. 

XLAM=XLAL 
RHO  =RHL 
XLAMB=XLA1 
XMUB-XM1 
RH0B=RH1 
WRITE  <  5 »  3 ) 

FORMAT (3X» 'LIQUID-SOLID  1  INTERFACE') 

WRITE (  5  »  4  )  THETAO 
FORMAT ( 3X » ' THETAO= ' » E12 . 5 ) 

CALL  WASO ( THETAO  >  XLAM » XMUB » XLAMB ,  RHO  f  RHOB ) 
CL=CLL 
CLB=CL1 

CALL  CRCSUB ( CL  *  THETAO » CLB  r  YB ) 

THET31=YB 
WRITE(5*5)  THET31 
FORMAT <3X» ' THETA31  =  ' *  El  2  *  5 ) 

LAYER  2 - SOLID  1-SOLID  2  INTERFACE 

XLAM^XLAl 
XMU-XMl 
RH0=RH1 
XLAMB=XLA2 
XMUB=XM2 
RHQB=RH2 
THETA0-THET31 
URITE«5r6) 

6  FORMAT <3X» 'SOLID  1-SOLID  2  INTERFACE') 

WRITE ( 5 »  7 )  THETAO 

7  FORMAT <3X» 'THETAO^' »E12.5) 

CALL  SOLID ( THETAO » XMU » XLAM > XMUi , XLAMB  r RHO  t RHOB ) 

CL=CL1 

CLB-CL2 
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CALL  CRCSUB ( CL »  7  HET  AO » CLB » YC  ) 

THET32-YC 
WRITE (5 » 8 )  THET32 
6  FORMAT ( 3X  »  ' THETA32= '  »E12»5) 

CALL  CRCSUB ( CL  »  THETAO » CL  >  YA  > 

THET12=YA 
WRITE (5*9)  THET12 
FORMAT <3X*  '  THETA  12=- '  *  E12 . 5 ) 

LAYER  3 - SOLID  1-LIGU1D  INTERFACE 

XLAM=XLA1 
XMU~XM1 
RHO-RH1 
XLAMB-XLAL 
RH05-RHL 
THET  AO-THET12 
WRITE  ('J»  10) 

F  ORMAT ■  3X  » ' SOLID  1-LIQUID  INTERFACE ' ) 

WRITE <5/ 11)  T  HE I AO 

11  FORMAT  <  3X  t  ' THET  AO  -  ' *  E12 . C  > 

CALL  SOLL I <  THETA 0 * XL AM r XML , XL  A  MB «  RHO  r RHOB ) 

C  LAYER  4 — SOLID  2-LIGU1D  INTERFACE 

XLAM-XLA2 
XMU-XM2 
RH0-RH2 
XLAMB-XLAL 
RHOB-RHL 
THETAO-ThEl 32 
WRITE ( 5  ? 12  > 

12  FORMAT ( 3X  * ’SOLID  2-lIQUID  INTERFACE' ) 
WR1TE(5»1S)  ! HET  AO 

13  FORMAT <3X» ' THETAO* ' * E12 . 5 ) 

CALL  SOLL  I  (  THET  AO *  XL  AM *  XMU  *  XLAiiB  f  RHO  t  RHOB ) 

CL  =  CL2 

CALL  CRCSUB ( CL  t THE TAG / CL * YA ) 

THET41=YA 

WRITE (5*14)  THET41 

14  FORMAT <3 X* ' THETA41-' >E12.5) 

C  LAYER  5 - SOLID  2-SOLID  1  INTERFACE 

XL.AM=XLA2 

XMU-XM2 

RHC-RH2 

XLAMB-XLwl 

XMUB-XMl 

RHOB-RHi 

THETA0=TMET41 

WRITERS* IS) 

15  FORMAT <  3X  * 'SOLID  2-SOLID  1  INTERFACE') 

WRITE (5* 16)  THET  AO 

16  FORMAT (3X* ' THETAO-  '*£12.5/ 

CALL  SOL  I D  ( THETAO  -  XhU »  XL  AM » XrlUfc  *  XuAKB  rhr.Qt  RHOB ) 

CL=CL2 

CL 6=0 LI 

-ALL  CRCSUB  (  CL  *  ■  i-.L )  i :  j  r  C;_;  *  i  C  / 

THE 735- YC 
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URI7E(5f 17)  THET35 

FORMAT  ( 3X »  '  THETA30-' '  r E12 . 5 ) 

LAYER  6 - SOLID  1-LIQUID  INTERFACE 

XLAM=XLA1 

XMU=XM1 

RHQ=RH1 

XLAMB=XLAL 

RH3B=RHL 

THETA0=THET30 

WRITE ( 5r 18 J 

FORMAT (3X» 'SOLID  1-LluUID  INTERFACE ' ) 

WRITE ( 5 » 19  y  in- I  AO 
F  OR HAT <  3X »  ' TKLT AO~  .£12.0) 

CALL  SOLL I  (  1  KET  AO  »  XuAr.  r  a:-iL  •  XLAMD  t  RHO  y  RHGB ) 

STOP 

END 
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SUBROUTINE  WASO ( THETAO  f XL AM  t XHUB  t XLAMB»  RHO t RHOB ) 

CL=SQRT ( XLAM*2 . 683/RHO ) 

CLB=SQRT( ( <XLAMB+2.*XMUB) ) *2 . 683/RHOB ) 

CTB=SQRT <  XMUB*2 . 683/RHOB ) 

CALL  CRCSUB<  CL  *  THETAO f  CL  f  YA ) 

THETA1=YA 

WRITE(5f2>  THETA1 

CALL  CRCSUB <  CL  f  THETAO  f CLB » YB  > 

THETA3=YB 

WRI TE ( 5  f  3 )  THETA3 

CALL  CRCSUB (CL  f  THETAO  f  CTB  f  YC ) 

TH£TA4=YC 

WRITE ( 5  f  4 )  THETA4 

2  FORMAT ( '  THETA1= ' » FIO . 4 ) 

3  FORMAT < '  TH£TA3=  '  f FIO . 4 ) 

4  FORMAT < '  TH£TA4= ' f  FIO ♦ 4 ) 

C11=-SIN( THETA 1 ) 

C12=SIN ( THETA3  > 

Cl 3= -COS ( THETA4 ) 

C21=C0S ( THETA 1 ) 

C22=C0S ( THETA3) 

C23=SIN(THETA4) 

C31=-XLAM/XLAM 

C32=CL* ( XLAMB+2 . *XMUB*COS < THETA3 ) *COS ( THETA3 ) ) / ( CLB*XLAM ) 
C33=(CL*XMUB*SIN(2.*THETA4) )/(XLAM*CTB) 

B1=SIN<  THETAO) 

B2=C0S(THETA0) 

B3=XLAM/XLAM 

WRITE  (5  f  1 )  Cll fC12fC13fB1 
WRITE (5f  1 )  C21  fC22fC23fB2 
WRITE ( 5  f 1 )  C31fC32fC33fB3 
1  FORMAT  <1XfE12«5f1XfE12*5f1XfE12»5f1XfE12.5) 

CALL  D3(DETERMfC11 fC12fC13fC21 fC22fC23fC31fC32fC33) 
WRITE(5r 10)  BETERM 
10  FORMAT ( 5X  f  ' DETERM= ' » E12 . 5 ) 

CALL  D3(BIfB1fC12fC13fB2fC22fC23fB3fC32fC33) 

X1=DI/BETERM 
WRITE(5f20)  XI 

20  FORMAT ( 5X  f ' A1/A0= ' f  E12 ♦ 5 ) 

CALL  D3(DII fC11fB1fC13fC21 fB2fC23fC31 fB3fC33) 

X2=DII/DETERM 

WRITE ( 5  f  30 )  X2 

30  FORMAT  <5Xf 'A2/A0=/ fE12*5) 

CALL  D3(DIIIfC11fC12fB1fC21fC22fB2fC31fC32fB3) 

X3=DI I I/DETERM 
WRITE<  5  f  40 )  X3 

40  FORMAT (5Xf 'A3/A0=' f£12»5) 

RETURN 

END 
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J. 


* 


-u  w  r  j 


SUBROUTINE  SOLLI ( THETAO f XLAM f XMU f XLAMB f RHO f RHOB ) 
CT=SQRT ( XMU*2 . 6B3/RH0 ) 

CL=SQRT  < ( <XLAM+2.*XMU)*2.683)/RHQ) 

CLB=SQRT ( XLAMB*2 . 683/RHO > 

CALL  CRCSUB(CLfTHETAOfCLfYA) 

THETA1=YA 
WRITE(5f2)  THETA1 
CALL  CRCSUB(CLfTHETAOfCTfYB) 

THETA2=YB 
URITE<5f3)  THETA2 
CALL  CRCSUB( CL f THETAO fCLBfYC) 

THETA3=YC 
WRIT£(5f4)  THETA3 
FORMAT ( '  THETAl=' fF10*4) 

FORMAT ( '  THETA2= ' f FI 0.4) 

FORMAT ('  THETA3= ' fF10«  4) 

Cl 1  =  -SIN( THETA1 ) 

C12=SIN(THETA2) 

C13=SIN( THE T A3) 

C21=C0S<  THETA 1 ) 

C22=-SIN( THETA2 ) 

C23=C0S(THETA3) 

B33=XLAM+2 .  *XMU*COS  (  THETAO  )  *CGS  <  THETAO ) 
C31=(XLAM+2.*XMU*C0S(THETA1 ) *COS < THETA1 ) )/E33 
C32=CL*XMU*SIN ( 2 . *THETA2 ) / ( CT*B33 ) 

C33=CL*XLAMB/ < CLB*B33 ) 

B1=SIN  < THET  AO ) 

B2=C0S( THETAO) 

B3= ( XLAM+2 . *XMU*COS ( THETAO ) *CQS (THETAO ) ) /B33 
WRITEOfI)  Cll  fC12fC13fB1 
WRITE (5f l )  C21fC22fC23fB2 
WRITE(SfI)  C31 fC32fC33fB3 
WRITE(5f5)  B33 
1  FORMAT (1XfE12«5f1XfE12*5f1XfE12.5f1XfE12«5> 

5  FORMAT  < '  B33= ' f  E12 ♦ 5 ) 

CALL  D3(DETERMfC11 fC12fC13fC21fC22fC23fC31fC32fC33) 
WRITE (5f 10)  DETERM 
10  FORMAT (5Xf 'DETERM®' fE12*5) 

CALL  D3< DI f B1 fC12fC13fB2fC22fC23fB3f C32 fC33 ) 

X1=DI/DETERM 

WRITE(5f20)  XI 

20  FORMAT ( 5X  f ' A1/A0= ' f  E12 ♦ 5 ) 

CALL  D3(DIIfC11fB1fC13fC21fB2fC23fC31 fB3fC33) 

X2=DII/DETERM 

WRITE(5f30)  X2 

30  F0RMAT(5Xf 'A2/A0=' fE12.5) 

CALL  D3(DIIIfC11fC12fB1fC21fC22fB2fC31fC32fB3) 

X3=DI I 1/DETERM 
WRITE ( 5  f  40 )  X3 

40  FORMAT (5Xf ' A3/A0® ' fE12*5) 

RETURN 

END 
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SUBROUTINE  SOL I D ( THETAO » XMU >  XLAM , XHUB  r  XLAHB » RHO  r  RHOB ) 
CL=SQRT ( (XLAM+2. *XMU)/RHO) 

CT=SQRT (XMU/RHQ) 

CLB=SQRT< ( XLAMB+2 . *XMUB  > /RHOB ) 

CTB=SGRT  ( XHUB/RHOB ) 

CALL  CRCSUB < CL » THETAO, CL fYA) 

THETA1=YA 

WRITE ( 5 » 1 )  THETA1 

CALL  CRCSUB (CL, THE T AO, CT  ,  YB) 

THETA2=YB 

WRITE <5, 2)  THETA2 

CALL  CRCSUB ( CL » THETAO , CLB » YC ) 

THETA3=YC 

WRI TE ( 5 » 3 )  THETA3 

CALL  CRCSUB ( CL , THETAO ,CTB »  YD ) 

THETA4=YD 

WRITE  (5,4)  THETA4 

1  FORMAT < '  THETA 1= ' » E12.5) 

2  FORMAT  < '  THETA2= ' , E12 . 5 ) 

3  FORMAT ( '  THETA3= ' , E12 . 5 ) 

4  FORMAT  ( '  THETA4= ' , E12 . 5 ) 

C11=-SIN(THETA1 ) 

C12=-C0S( THETA2) 

C13=SIN(THETA3) 

C14=-C0S(THETA4) 

C21=C0S(THETA1) 

C22=-SIN ( THETA2 ) 

C23=COS (THETA3) 

C24=SIN(THETA4) 

B33=XMU 

C31=XMU*SIN<2.*THETA1)/B33 

C32= ( XMU*CL*COS ( 2 . *THETA2 ))/( C T*B33 ) 

C33=XMUB*CL*SIN  < 2 . *THETA3 ) / < CLB*B33 ) 
C34=-XMUB*CL*COS(2.*THETA4)/(CTB*B33) 

B44=XLAM+2 . *XMU*CGS  <  THETAO  >  *COS ( THE 7  AO ) 

C41=- ( XLAM+2 . *XMU*COS ( THETA 1 ) *COS < THETAl ) > /B44 
C42=XMU*CL*SIN(2,*THETA2)/(CT*B44) 

C43=CL*( XLAMB+2. *XMUB*COS < THETA3 ) *COS ( THETA3 ) )/(CLB*B44) 
C44=XMUB*CL*SIN(2.*THETA4)/(CTB*B44> 

B1=SIN ( THETAO ) 

B2=CQS( THETAO) 

B3=XMU*SIN(2. ♦THETAO ) /B33 

B4= ( XLAM+2 . *XMU*COS  <  THETAO ) *COS  <  THETAO ) ) /B44 

WRITE (5*5)  C11,C12,C13,C14,B1 

WRITE ( 5 » 5 )  C21,C22,C23,C24,B2 

WRITE (5 ,5)  C31 , C32 , C33 » C34 , B3 

WRITE ( 5 , 5 )  C41,C42,C43»C44,B4 

5  F0RMAT(1X,E12.5,1X,E12.5,1X,E12.5,1X,E12.5, 1X,E12.5) 

CALL  D4  <  DETERM » Cl 1  ,C12,C13,C14,C21,C22,C23,C24»C.31, 

XC32 , C33 , C34  »C41»C42,C43,C44) 

WRITE <5, 10)  DETERM 

10  FORMAT  ( 5X »  ' DETERM= ' , El  2 . 5 ) 
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CALL  D4<DI,BltC12fC13rC14»B2»C22fC23»C24, 

XB3 1 C32 1 C33  r  C34  fB4»C42»C43  *C44 ) 

X1=DI/DETERM 
WRITE (5»20>  XI 

20  FORMAT  <5Xf /A1/A0=/ f E12*5) 

CALL  D4(DII»Cll»Bl»C13fC14,C21»B2rC23rC24»C31f 
XB3 » C33 1 C34  fC41»B4»C43*C44) 

X2=DI I/DETERM 
WRITE < 5 »  30 )  X2 

30  FORMAT  <  5X  f  ' A2/A0= ' » E12  *5 ) 

CALL  D4<DIIIfC11fC12fB1fC14fC21fC22fB2fC24f 
XC31 fC32fB3fC34fC41 fC42fB4fC44> 

X3=DI I I/BETERM 
WRITE<5f40)  X3 

40  FORMAT ( 5X f  ' A3/A0- ' f  E12 . 5 ) 

CALL  D4(DIVfC11fC12fC13fB1fC21fC22fC23fB2fC31f 
XC32fC33fB3fC41fC42fC43fB4) 

X4=DIV/DETERM 
WRITE<5f50)  X4 

50  FORMAT (5X  f ' A4/A0= ' »E12.5) 

RETURN 

END 
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SUBROUTINE  D4 ( D > A1 1 » A12 » A13 » A1 4 , A21 > A22 , A23 t A24 » h31 t A32 r 
XA33»A34»A41*A42»A43>A44) 

CALL  D3 ( DA » A22 >A23>A24»A32»A33» A34 »A42*A43»A44) 

CALL  D3(DB> A21 > A23 > A24 » A31 » A33 > A34 r A41fA43>A44) 

CALL  D3<DC,A21,A22>A24,A31» A32,A34,A41> A42»A44> 

CALL  D3 ( DD » A21 » A22 1 A23 »A31»A32»A33»A41jA42>A43) 
D=A11*DA-A12*DB+A13*DC-A14*DD 
RETURN 
END 
C 

SUBROUTINE  D3(D»A11»A12»A13»A21»A22>A23»A31 , 

XA32  r A33 ) 

CALL  D2<DE»A22r A23»A32»A33) 

CALL  D2(DF»A21 ,A23, A31 »A33> 

CALL  D2 ( DG »A21»A22»A31»A32> 

D=A11*DE-A12*DF+A13*DG 

RETURN 

END 

C 

SUBROUTINE  D2 < D t A1 1 r A12f A21 t A22) 

D=A11*A22-A12*A21 

RETURN 

END 

* 


SUBROUTINE  CRCSUB < A , B » C t Y ) 

X=A*SIN(B)/C 

Y=X+(X**3. >/6.+3.*<X**5. )/40,+5.*(X**7. )/112. 
X+33.*<X**9. )/576. 

RETURN 

END 
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APPENDIX  B 

TWO-DIMENSIONAL  CORRELATION 
OF  TWO  LINES  X  =  Y  =  UNIFORM  TRANSLATION 
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Li  I  MENS  I  UN  i  J  00  »  h.i  ■  i  A  a  v  j  O )  r  m  Y  <  J  O  /  •  til  1 1 1  i  00  ) 
RE  A  Li  (  7 » 1 0  >  H  *  i'i  >  Nf 
READ  (  / »  1 1  >  L'ELX 
:  0  F OPhA K w .  .  / 

j.  1  F  QRhA  I  .  F  i  kj  ,  4  ) 

Fill  1  I  j  «  N  » i 
DO  1.3 

...  A  < 1 1 J  >  •  0 »  c 

3,:  CONTINUE. 

READ  v  7,20;  ■.  D I  (  i  »  w  1 1  I » I 00  ) 

.0  F ORKAT ( t  10 . A ) 

DO  3  I - 1 » 1 uO  »  L 

A  ( 1 » 1 )  -  D  Kill) 

3  CONTINUE 

DO  5  1=1.99,1 
DO  6  J=1 , 99  » 1 

6  A <1  +  1 . J! 1 >  -A< 1  +  1  *  1 )*A(lf J+l)/3. 

5  CONTINUE 

N=N~M 

DU  9  Mti-0  » K,  M 
NN=MM 

DO  7  1=1 f M» 1 

AX  (  I )  ••  •  A  <  MM  + 1  f  I +NN  ) 

AY ( I ; =A  < I  +  MM »  NN+1 ) 

7  CONTINUE 

CALL  CORK' I <  A  r  AX , A  Y , N . M , NP  > 

9  CONTINUE 

STOP 
END 
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c 

SUBROUT IN* h;  ClIR. .  i  (  A  t  AX  r  M  l  r  N  r  H  f  IC  ) 

WR 1 1  E  ( b * 30 > 

30  F  ORMA  i  (  OX  r  '  COREl  > 

D  I M  I::.  N  S  1  u  N  A  i  i  0  0  >■  i.  00  t  >  A  X<  10)  t  A  Y  ( 1 0 ) 

S2-- 1  0E20 

C2~  I0i£20 

INM-:0 

!.  H  N  ■  ■  0 

IwNM-N-M 

L‘3  2  J~i  t  ICNi'i  *  1 

DO  1  L  •  0 >  ICNM » 1 

S-0.0 

C --  0  .  0 

DO  3  1 1-1 »  M  » i 
SI -AX  <  II )  ~A  <.J » 1 1  tl.  > 

C  .1  -AY  \  1 1 ) -A  < 1 1 +L  r  J ) 

S-=S  +  S1*S1 
C=C+C1*C1 
3  CONTINUE 

IF" <  S .  LT .  S2  >  INM--L 
IF (S»LT *S2)  S2--S 
IF  t  C . LT .  C2  <  IMN-L 
IF (C*LT *C2)  C2-C 

1  CONTINUE 

2  CONTINUE 

DX=0 *001 *FLOAT  < I NM* I C ) 

DY^O*  001*  FLOAT  <IMN*IC) 

WRTTE<5»?)  S2 

9  FORMAT <5X» ' CORRELATION-X- ' »E12*5) 

WRITE ( 5 » 7 )  C2 

7  FORMAT <5X* ' CORREALTION-Y= ' r E12 . 5 ) 
WRITE < 5* 8 )  DX > BY 

8  FORMAT  <:  OX »  ■' DX  - '  » FI 0 . 3  >  3X ,  '  DY= '  » F 10 . 3 ) 
RETURN 

END 
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CORRI 

CORLELnaON-X- 
CORREAL! 10N-Y- 

nx~  o.ooo 

CO  RFC! 

CORRELA! 1GN~X~ 
CORREAL! ION-Y= 
DX®  O.005 

GORKI 

COR R £ !...  AT  1 0 ':•> -  X  - 
CORREA...  !  CON- Y- 
I»X=  0.0  iO 
CORRI 

CORRELA  T I ON -X- 
CuRIxFJ AL  »  iON-Y« 
DX-  0.0 IS 

CURRI 

C 0 RREL A  I  1 0 N  -X~ 
CORREAL! I ON- Y= 
DX=  0.020 

CORRI 

COR RELAY IQN~X~ 
CORREAL  I 10N-Y® 
DX----  0.025 

CORRI 

CORRELA  TIGN-X® 
CORREAL !  YLGN-Y- 
DX=  0.030 

CORRI 

CORRELA! ION- X= 
CORREALTION-Y® 
DX=  0.035 

CORRI 

CORRELA rroN-x= 

CORREALTION-Y- 
1)X~  0.040 

CORRI 

C  0  R  R  E  L  A  f  x  0  iV — X  — ' 
CORREALTION-Y= 
DX-  0,045 

CORRI 

CORRELA TION-X= 
CORREAL! ION-Y- 
DX-  0.050 


0 .  vOvvvt-.r  Ou 
0 . 0 0000c TOO 
DY -  O.uOO 

0 . OOOOOE+OO 
0 » OOOOOE+OO 
{(:-•-  0,005 

0 » OOOOOE+OO 
0 .  OOOOOE+OO 
i‘Y-  O.OlO 

0. OOOOOE+OO 
0 . OOOOOE+OO 
DY-  0.015 

0 . OOOOOE+OO 
0 . OOOOOE+OO 
D  Y  ~  0.020 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY-  0.025 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY-  0.030 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY=  0,035 

0, OOOOOE+OO 
0. OOOOOE+OO 
DY=  0.040 

0  - OOOOOE+OO 
0. OOOOOE+OO 
DY-  0.045 

0, OOOOOE+OO 
0. OOOOOE+OO 
D i  0,050 
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CORF:] 

CQRRELA  I  ION-X= 
CORREAL [ ION  - Y~ 
D X  =  0.055 

CORRI 

CORRELATION— X= 
CORREALTION-Y* 
DX  =  0.060 

CORRI 

CORRELATION- X= 
CORREALTION-Y - 
DX~  0.065 

CORRI 

CORRELATION- X- 
CORREALTION-Y* 
DX=  0.070 

CORRI 

C0RRELA1 ION-X- 
CORREAL.  I ION-Y- 
»X*  0.075 

CORRI 

CORRELATIONS* 
CORREALTION-Y* 
DX=  0.080 

CORRI 

CORRELATION-X= 
CORREALTION-Y* 
DX*  0.085 

CORRI 

CORRELATION- X* 
CORREALTION-Y* 
DX  =  0,090 

CORRI 

CORRELATION-X- 
CORREALTION-Y- 
DX-  0.095 


0  •  +00 
0.0 OOOOE+OO 
LA-  0.055 

0 . OOOOOE+OO 
0.05000E+00 
HY  •  0.060 

0 . OOOOOE+OO 
0 , OOOOOE+OO 
IiY-  0.065 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY*  0.070 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY*  0.075 

0. OOOOOE+OO 
0. OOOOOE+OO 
Li  Y  *  0.080 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY*  0.085 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY-  0.090 

0. OOOOOE+OO 
0. OOOOOE+OO 
DY-  0.095 
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